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ABSTRACT 

We report new spectroscopic and photometric observations of the parent stars of 
the recently discovered transiting planets TrES-3 and TrES-4. A detailed abundance 
analysis based on high-resolution spectra yields [Fe/H] = —0.19 ± 0.08, T^d = 5650 ± 
75 K, and log g = 4.4 ± 0.1 for TrES-3, and [Fe/H] = +0.14 ± 0.09, T^s = 6200 ± 75 K, 
and log 5 = 4.0ib0.1 for TrES-4. The accuracy of the effective temperatures is supported 
by a number of independent consistency checks. The spectroscopic orbital solution for 
TrES-3 is improved with our new radial- velocity measurements of that system, as are the 
light-curve parameters for both systems based on newly acquired photometry for TrES-3 
and a reanalysis of existing photometry for TrES-4. We have redetermined the stellar 
parameters taking advantage of the strong constraint provided by the light curves in the 
form of the normalized separation a/ii^ (related to the stellar density) in conjunction 



^Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138 USA; 
asozzett@cfa.harvard.edu 

^INAF - Osservatorio Astronomico di Torino, 10025 Pino Torinese, Italy 

■^Department of Physics, and Kavh Institute for Astrophysics and Space Research, Massachusetts Institute of 
Technology, Cambridge, MA 02139 USA 

^Department of Physics & Astronomy, Bowling Green State University, Bowling Green, OH 43403 USA 

^California Institute of Technology, 1200 East California Boulevard, Pasadena, CA 91125, USA 

''NASA Postdoctoral Program Fellow, Goddard Space Flight Center, 8800 Greenbeh Rd Code 690.3, Greenbelt 
MD 20771 USA 

'^Lowell Observatory, Flagstaff, AZ, USA 

^Planetary Science Institute, Tucson, AZ 85719 USA 

^Instituto de Astrofi'sica de Canarias, 38200 La Laguna, Tenerife, Spain 
^"Las Cumbres Observatory Global Telescope, Goleta, CA 93117 
^^Department of Physics, University of California, Santa Barbara, CA 93106 
Alfred P. Sloan Research Fellow 



- 2 - 



with our new temperatures and metallicities. The masses and radii we derive are 
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for TrES-3 and TrES-4, respectively. With these revised stellar parameters we obtain 
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for TrES-4. We confirm TrES-4 as the planet with the largest radius among the currently 
known transiting hot Jupiters. 



improved values for the planetary masses and radii. We find Mp 
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Subject headings: stars: individual (TrES-3) — stars: individual (TrES-4) — stars: 
abundances — stars: fundamental parameters — planetary systems 



1. Introduction 

The 40 transiting planet systems confirmed as of August 200^ show a remarkable diver- 
sity of properties, which is indicative of the complexity of planet formation and evolution pro- 
cesses. Many different follow-up studies enab l ed by the special orientation of these systems ( e.g. 



Queloz et al 



200C 



Charbonneau et al. 



2001 



2005 



Knutson et al 



2007 



Tinetti et al. 



20071 . see 



Charbonneau et al. 2007 for a review) have br ought about rapid improvements in evolution- 



ary models of planet interiors and atm ospheres ([Baraffe et al 



2004 



Lecavelier des Etangs 



2006 



Guillot et al.ll2006l : iBurrows et al.ll2007l ). The increasing predictive power of these models is begin- 
ning to drive even more challenging observations. On the other hand, the precision and accuracy 
with which the most basic planet properties such as the mass and radius can be determined is 
currently limited by our knowledge of the properties of the host stars. Significant uncertainties 
remain in the stellar mass and radius determinations of many systems. In some cases this is due to 
poorly determined photospheric properties (mainly temperature and metallicity) , and in others to 
a lack of an accurate luminosity estimate. Additionally, the variety of methodologies used for these 
determinations and the different approaches towards syst ematic errors have resulted in a rather 
inhomogeneous set of planet properties, as discussed by [Torres et al.l (|2008l ). This complicates 
the interpretation of patterns and correlat ions that are being proposed (e.g., iMazeh et al.ll2005l : 
Guillot et al.ll2006l : [Hansen &: Barmanll2007l ). Recent improvements in the analysis techniques have 
the potential to increase the accuracy of the stellar and planetary parameters significantly, espe- 
cially for the large majority without a direct distance estimate. In particular, the application of the 
constraint on the stellar density that comes directly from the transit light curves has been shown to 
be superior to the use of other indicato r s of luminosity such as the surface g r avity (log 5) determined 



spectroscopically (jSozzetti et al 



2007 



Holman et al.l 120071 ). iTorres et al.l ([2003) have recently re- 



analyzed a large subset of the known transiting planets, incorporating these improvements and 
applying a uniform methodology to all systems. 



^For a comploto listing seej http: //www. inscience . ch/transitsT] or |http : / /exoplanet . eu] . 
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In the present wor k we focus on two of the recently dis covered transiting systems, TrES-3 
(jO'Donovan et alj 120071 ) and TrES-4 (jMandushev et al.l 120071 ) . which lack accurate estimates for 
the photospheric properties of the parent stars and as a result have more uncertain stellar and plane- 
tary parameters. To improve upon these properties, we present new radial- velocity and photometric 
observations of TrES-3 with which we refine both the light curve solution and the spectroscopic 
orbit. We also carry out a reanalysis of existing TrES-4 photometry utilizing a technique which 
treats stellar limb darkening using adjustable parameters. We perform the first detailed spectro- 
scopic determination of the photospheric properties of both stars, and we make use of the constraint 
on the stellar density mentioned above to infer more accurate values for the stellar and planetary 
masses and radii. Our paper is organized as follows. In §[2] we summarize the observations. In §[3] 
we present our effective temperature and metallicity determinations, along with several consistency 
checks aimed at establishing the accuracy of the temperatures. In §2] we report an updated spec- 
troscopic orbital solution for TrES-3, and the light curve solutions for both systems are discussed 
in §[5l Section[6] then describes our determination of the stellar masses and radii, which in turn 
lead to refined values for the planetary parameters over those reported in the discovery papers. We 
conclude in §[7] by providing a summary of our results and by discussing whether the properties of 
the host stars give any useful clues on the origin of the strongly contrasting densities of their close-in 
gas giant planets, particularly in comparison with the other known transiting planet systems. 



2. Observations 



2.1. New radial velocities for TrES-3 



High-resolution, high-SNR (signal-to-noise ratio) spectroscopic observations of TrES-3 were 
obtained in July and September 2007 with the HIRES spectrograph on the Keck I telescope 
(jVogt et al.lll994l ). using essentially the same setup as in the discovery paper. The spectra cover 
the effective wavelength range ~3200-8800 A, and an iodine (I2) gas absorption cell placed in front 
of the spectrograph slit was used to superimpose a dense set of narrow molecular lines between 
~5000 A and ~6000 A. This iodine spectrum provides a stable wavelength reference and a means 
of monitoring the instrumental profile that are crucial for achieving precise radial velocity determi- 
nations. The observations were reduced and extracted using the MAKEE software written by T. 
Barlow. We obtained 4 spectra with typical exposure times of 15 min re sulting in SNRs o f ~100 



pixel ^, which we combined with the discovery observations reported by lO 'Donovan et al.l ( 20071 ) 
for a total of 11 iodine spectra. 

Precise radial velocities were measured following a p rocedure based on the methodology devel- 
oped for the AFOE spectrograph (jKorzennik et al.ll200d ). and adapted for the processing of HIRES 
spectra. Conceptually, the method models the observed star and I2 spectrum by using templates 
for the stellar and I2 spectra. This model includes physically motivated parameters that describe 
the spatial and temporal variations of the instrument, including the instrumental profile, as well as 
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the sought-after relative Doppler shift with respect to the stellar template. The model parameters 
are adjusted to minimize the difference between the model and the observations in a least-squares 
sense, down to the Poisson noise. The stellar template is estimated by deconvolving an observation 
taken without superimposing I2, whereas the template for the I2 is based on a high-resolution and 
high SNR scan of the gas cell using the F ourier Transform S pectrometer on the McMath Solar 
Telescope at the Kitt Peak Observatory (see lButler et al.lll996l ). Each spectral order is modeled as 
a whole but analyzed independently. The final radial velocities are estimated from the mean of the 
Doppler shifts computed for each order and their uncertainties from the standard deviation of that 
mean. The resulting velocities, expressed in the solar system barycentric frame, along with their 
associated formal errors are reported in Table [TJ These measurements include and supersede the 
velocities presented in the discovery paper, and they also correct a minor error in the previously 
published dates of observation. The revised orbital solution is presented and discussed in §|H 



2.2. Differential photometry: new TrES-3 data and revisited TrES-4 data 



In addition to the photometric measurements presented in the discovery paper (lO 'Donovan et al 



20071 ). we have collected six other high-cadence, high-precision transit light-curves of TrES-3. One 
was obtained in the V band on UT 2007 April 23 using the CCD camera of the IAC80 telescope 
at the Ob servatorio del Teide , Tenerife, Spain, and the other five were gathered using Kepler Cam 



(see, e.g., iHolman et al.ll2006l ) at the FLWO 1.2-m telescope: two in the Sloan g and r bands on 
UT 2007 April 25 and 2008 March 27, and the remaining three in the i band on UT 2008 March 
9, April 12, and May 8. On all nights conditions were quite good, except for the presence of high 
cirrus and highly variable seeing on UT 2008 April 12. The transits were observed at airmasses 
ranging between ~1.00 and ~1.80. 

All datasets were reduced using standard calibration techniques (overscan correction, trimming, 
bias subtraction, flat fielding). We then performed aperture photometry of TrES-3 and between 
10 and 30 comparison stars, depending on the filter and exposure time. We experimented with 
different choices for the aperture size, comparison star ensemble, and weighting of the comparison 
stars, aiming for the smallest out-of-transit (OOT) RMS. In practice, the best aperture size was 
approximately twice the FWHM of the stellar image on each night, and the best results were 
obtained using a straight average of the normalized light curves of the comparison stars. As for the 
formal errors on each photometric data point, we used the product of the OOT RMS of each light- 
curve and the f actor > 1, used t o account for departures from Gaussian ( "white" ) uncorrelated 



noise (see, e.g.. lWinn et al.l l2008bl ) . 



For consistency, we also re-processed in the same way the two discovery light-curves of TrES-3 
presented in O'Donovan et al. (2007), as well as the three photometric datasets utilized by Man- 
dushev et al. (2007) in their discovery announcement of TrES-4. 

The main characteristics of the TrES-3 and TrES-4 light-curves are summarized in Table [2j 
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The final set of photometric time-series of TrES-3 in all filters (including the discovery data) is 
available in machine-readable form in the electronic version of Table [3l and is plotted in Figure [TJ 
Table U] reports all photometric data for TrES-4, which are shown graphically in Figure [2j 



3. Atmospheric parameters and age constraints 

3.1. Spectroscopic abundance analysis 

A detailed abundance analysis was carried out using the Keck/HIRES template spectra of 
TrES-3 and TrES-4. The stellar a tmosphere parameters Tpff, log q, and [Fe/H] were determined 



using standard methodology (e.g., iGonzalez &: LambertI Il996l : iGonzalez et al.l l200ll : I Santos et al 



20041 ). which we summarize as follows. Equivalent widths (EWs) for a set of relatively weak Fe I 
and Fe II lines were mea sured manually in the Keck spectra using the splot task in IRAF (see, 



e.g.. iSozzetti et al.ll2004l . and references therein, for details on the specific choice of lines). These 
were then use d together with a grid of Kurucz ATLAS plane-parallel stellar mo del atmospher es 
( Kurucz 1993 ) as inputs to the 2002 version of the MOOG spectral synthesis code ( Sneden 1973l fl 
Atmospheric parameters were derived under the assumption of local thermodynamic equilibrium 
(LTE), imposing excitation and ionization balance. Formal uncert ainties on Tp^, logq, and in icro- 
turbulent velocity Et were d erived using the approach described in iNeuforge MagainI (jl997l ) and 
Gonzalez &: Vanturd (jl998l ). while the nominal uncertainty for [Fe/H] corresponds to the scatter 
obtained from the Fe I lines rather than the formal error of the mean. The resulting set of pa- 
rameters for TrES-3 is T^s = 5650 ± 75 K, log 5 = 4.4 ± 0.1, ^ = 0.85 ± 0.05 kms-\ and [Fe/H] 
= -0.19 ± 0.08, and for TrES-4 we obtain T^^ = 6200 ± 75 K, logg = 4.0 ± 0.1, Et = 1-50 ± 0.05 
kms-\ and [Fe/H] = -F0.14±0.09. These values are collected in Tableland Table [HI respectively. 

As a check, the manual EW measurements used above were com pared against the re sults from 
the automated software AREi , made available to the community by ISousa et al.l (120071 ) . Figure [3] 
shows the comparison between the EWs measured for TrES-3 using ARES and those from the 
manual approach with IRAF. The top panel indicates excellent agreement between the two; the 
linear fit has a slope of 1.011. There are no significant correlations with the Fe I line strength 
(see middle panel): the fractional difference between EWs measured by the two methods has a 
probability of no correlation of 0.1. Finally, the histogram in the bottom panel shows there is no 
appreciable systematic difference between the methods. The mean diff erence and s c atter are only 
0.67 mA and 2.6 mA, respectively, which are smaller than found by ISousa et al.l (|2007l ) from a 
similar comparison between FEROS, HARPS and UVES spectra. Consequently, the results of the 
abundance analysis with MOOG using EWs measured with ARES are virtually identical to those 
presented earlier. 



^http : //verdi . as .utexas . edu/moog.htinl| . 
• ^http: //www, astro .up.pt /~ sousasag/ares| . 
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Following iGonzalea (j2008l ) we also synthesized a number of unblended Fe I lines in the template 
spectra of both TrES-3 and TrES-4, and determined projected rotational velocities of vsini = 
1.5 =b 1.0 kms~^ and vs'mi = 8.5 =b 0.5 kms~^, respectively. 

Our new determinations of the atmospheric parameters for th e parent stars of IVES-S and 
TrES-4 are genera l ly in good agreement with the values reported by lO 'Donovan et al.l (120071 ) and 
Mandushev et al.l (|2007l ). the main difference being that those authors assumed [Fe/H] = 0.0 in 
their studies, whereas we find significant departures from solar metallicity in both stars. 



3.2. External checks on T^g 

We describe here additional estimates of the effective temperature for TrES-3 and TrES-4 that 
serve to test the accuracy of our determinations above. 



3.2.1. CfA spectroscopy 

Both planet host stars were obs erved spectroscopically with the Center for Astrophysics (CfA) 
Digital Speedometer (jLathamlll992l ) as part of the regular follow-up after the discovery of photo- 
metric signals suggesting transits. These observations cover 45 A in a single echelle order centered 
at 5187 A, and have A/AA 35,000. While this resolving power is moderately high, the SNRs of 
the spectra are low, ranging from 7 to 13 per resolution element of 8.5 kms~^ for the 13 exposures 
of TrES-3, and SNRs of 11-13 for the 7 exposures of TrES-4. Nevertheless, useful information 



on the stellar properties can be extracted from these spectra as described by lTorres et al.l (|2002l ). 
Briefiy, the observed spectra a re cross-correlated against a library of syn thetic spectra based on 
Kurucz model atmospheres (see lNordstrom et al.lll994l : iLatham et al.ll2002l ). calculated over a wide 
range of values of Teg, log 5, [Fe/H], and rotational velocity Vrot- The combination of parameters 
yielding the highest correlation averaged over all exposures is adopted as a representation of the 
properties of the star. Due to the narrow wavelength range and limited SNRs, the first three of 
the above properties are typically strongly correlated and are difficult to determine simultaneously. 
Given these constraints, in the discovery papers we initially held the metallicity fixed at the solar 
value for both stars. 

For TrES-3, which our Keck spectroscopy indicates is slightly metal-poor, we repeated the 
above determination using a fixed metallicity of [Fe/H] = —0.5, and interpolated T^s and log 5 to 
the precise composition from Table[71 We obtained Tgff = 5530ibl30 K and log g = 4.45ib0.17, where 
the uncertainties include contributions both from the error in the adopted [Fe/H] and the internal 
errors. The temperature as well as the surface gravity are consistent with the Keck determinations, 
within the uncertainties. Similarly for TrES-4 we repeated the determinations for [Fe/H] = -1-0.5, 
given the metal-rich composition indicated by the Keck spectroscopy, and interpolated to the 
intermediate value from TableO The results are Tgff = 6270 it 150 K and logg = 3.96 it 0.17, which 
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are once again in good agreement with the more reliable estimate from Keck. 



3.2.2. Line depth ratios 



While the temp erature determinat i ons d iscussed so far are implicitly based on the strength 
of the spectral lines, iGray &: JohansonI (jl99ll ) have demonstrated that highly precise information 
can also be extracted using the ratio of the depths of two spectral lines having different sensitivity 
to temperature. However, rather than yielding absolute temperatures, this technique in its origi- 
nal formulation only measures c hanges in temperature, albeit with extremely high prec ision often 



reaching a few Kelvin (see, e.g.. lGravlll994l : ICatalano et al. 



2002 



Kovtvukh et al.ll2003l ^. Absolute 



temperatures can still be obtained with recourse to external color-temperature calibrations, since 
line-depth ratios (LDRs) are usually strongly correlated with the color index of the star. The 
accuracy of such T^s determinations is then limited by the calibrations themselves. 



Biazzo et al.l (|2007l ) have presented LDR-Tefj calibrations based on 26 carefully selected lines 



of Fe, V, Sc, Si, and Ni in the spectral interval 6190-6280 A, grouped into 16 line pairs. These cali- 
brations are valid for stars with temperatures between ~3800 K and ~6000 K, so unfortunately the 
technique is not applicable to TrES- 4. For TrES-3 we in easured all 26 lines in our high-SNR Keck 
template spectrum, and adopted the lBiazzo et al.l (j2007l ) calibrations app ropriate for stars with ro - 
tational velocities of 10 kms~^ (very close to the vsini value reported bv lMandushev et al.l ( 20071 ) 
and obtained also in this wo rk) . The tem perature scale of these relations relies on a transformation 
between B — V and Tgfj by iGravl (|2005l ) , which makes use of a mixture of dwarf and giant tem- 
peratures obtained by many different methods and does not account for differences in metallicity. 
For the pre sent work we have pre f erred to us e more sophisticat e d colo r-temperature relations such 
as those by lRamirez Melendea (j2005l ) and ICasagrande et al.l (|2006l ). which account not only for 
luminosity class but include also metallicity terms, and are based on effective temperatures derived 
homogen eously by the Infra red Flux Method. The conversion from the LDR-based T^ff i n ferred 
from the iBiazzo et al.l (j2007l ) relations back to an average color for the star using the iGravl (|2005l ) 
prescription gives B—V = 0.641 ±0.007. Our two preferred color-temperature relations mentioned 
above then yield a weighted average temperature of T^fj = 5710 it 70 K for TrES-3, in which the 
uncertainty includes observational errors propagated from the LDR measurements as well as the 
scatter of the calibrations. This result is consistent with our more direct estimate in § 13.11 



3.2.3. Ha line profiles 



As is well known, the wings of the line (but not its core, formed higher up in the atmosphere 
under non-LTE conditions) are ver y sensitive to Tpff variations, but are relatively insensitive to 
changes in log g and [Fe/H] (see, e.g.. lSozzetti et al. ISantos et ahlboO^ . and references therein) . 
This allows for a useful consistency check on our T^s estimates above. We compared the observed 
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Ha line profiles in our Keck template spectra of TrES-3 and TrES-4 against synthetic profiles for 
solar-metallicity dwarfs ([Fe/H] = 0.0, log 5 = 4.5) from the Kurucz database. The results of this 
exercise are displayed in Figure [U in which 10 A regions centered on Hq, are shown for each star 
together with four calculated profiles for different values of Tefj . In both cases the temperatures one 
would infer from these comparisons agree well with the estimates reported in Table [7| and Table [8j 



3.2.4- Photometric estimates 



An additional check on the effective temperatures is available from the multi-color photom- 
etry for TrES-3 and TrES-4. Measurements in Johnson BV, Cousins RI, and 2MASS JHKg 
were used to derive seven color indices for TrES-3, and nine for TrES-4 when considering also 
the Bt and Vr measureme nts fr om the Tycho-2 catalog (jH0g et al.l I2OOOI ) . The calibrations by 
Ramirez h Melendea (j2005l ) and ICasagrande et al.l (|2006l ) then yielded weighted average temper- 
atures of 5390 lb 50 K and 6135 it 50 K for TrES-3 and TrES-4, respectively, ignoring extinction 
and adopting [Fe/H] in each case as determined in §[3Tj The uncertainties include photometric 
errors and metallicity errors, as well as the scatter of the calibrations, but exclude unquantified 
systematic errors in the calibrations themselves. These temperatures are 260 K and 65 K cooler 
than our spectroscopic determinations above. They can be reconciled with the values in § 13.11 if we 
assume the presence of reddening, and correct each of the indices. For TrES-4 the required value of 
E{B-V) is hardly significant (0.013±0.010 ma g), but for TrES-3 w e obta in E( B-V) = 0.071 i fc0.013 
mag0 For comparison, the reddening maps of Burstein &: Heiles ( 1982 ) and Schlegel et al. ( 1998 ) 
indicate a total reddening along the line of sight to each star of only E{B — V) ^ 0.028 mag, and 



simila r results are obtained from the model of Galactic dust distribution by iDrimmel &: Spergel 
(j200ll ). This is not inconsistent with the small value we infer for TrES-4, but it is much smaller 
than our estimate for TrES-3. Possible explanations include patchy inters tellar material comb ined 
with the relatively coarse resolution of these maps (a few arc minutes for lSchlegel et al.lll998l ). or 
perhaps the presence of circumstellar material in TrES-3. 



3.3. Constraints on the age 



The reliability of age indicators for stars older than 1-2 Gyr, such as chromospheric activity and 
lithium (Li) ab undance, as well as their interagreement, have been the subject of much debate in the 



literature (e.g., Pace Pasauini 2004 ; Lambert fc Reddvir2004 ; Song et al. 2004 ; Sestito Randich 



2005; 



Sestito et al 



20061 ). The difficulties are due to many factors, including (but not limited to) 



non-trivial correlations between chromospheric activity, rotation, mass, and age, limited availability 



^The presence of reddening in TrES-3 is already apparent from our results in iil3.2.2l in which the LD R-based 
B~V color is significantly bluer than the measured value of B~V — 0.712 ± 0.009 IjO'Donovan et al.ll2007r ). 
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of activity estimates averaged over entire stellar activity cycleJl, insufficient understanding of the 
temporal evolution of Li depletion due to the complex interplay between various processes (e.g., 
convection, mixing, diffusion, mass loss), and non-negligible differences in the observed behavior 
of chromospheric activity and Li depletion as a function of mass, age, and chemical abundance 
between stars in young and old clusters, and in the field. 

We have nonetheless attempted to use the Ca II activity indicator and the lithium abundance 
as measured in our HIRES spectra of TrES-3 and TrES-4 to provide independent constraints for 
comparison with the formal age estimates determined below from evolutionary models (see §[6]), as 
well as to search for possible peculiarities of these two planet hosts compared to other samples of 
stars with and without planets. 

Figure [5] shows a region of the HIRES template spectra for each star centered on the Ca II 
H line. Clear emissio n is seen in th e core of the line of TrES-3, but not in TrES-4. Following the 
procedure outlined in lSozzetti et al.l (|2004l ) we measured the chromospheric emission ratio log R'fjx^ 
corrected for the photospheric contribution, from the Ca II H and K lines in each of our spectra. 
We obtained (logi?^^) = -4.54 ± 0.13 and (logi?^^) = -5.11 ± 0.15 for TrES-3 and TrES-4, 
respectively. TrES-3 thus appears moderately active, while TrES-4 is quite the opposite. The 
resulting chromospheric age estimates, using the relations summarized in lwright et al.l ( 2004 ). are 
t = 0.9 ± 0.7 Gyr for TrES-3 and t = 9.4 ± 1.7 Gyr for TrES-4. 

In Figure [6] we show the results of the spectral synthesis of a 10 A region centered on the 
Li A67 '07.8 line, using the atmospheric parameters derived from the Fe-line analysis and the line 
list of iReddv et al.l (|2002l ). The two panels display the comparison between the observed spectra 
and three synthetic spectra, each differing only in the assumed Li abundance. Neither star shows 
a measurable Li line, and we can only place upper limits of log e (Li) < 1.0 and log e (Li) < 1.5 for 
TrES-3 and TrES-4, respectively. By comparison with average Li abundan ce curves as a function of 
effective temperature for clusters of different ages (|Sestito fc Eandichlbood ). one would infer a rather 
old age for TrES-3 of t > 4 Gyr, nominally inconsistent with the estimate above from the Ca II 
activity index. Based on its Li the object a ppears decidedly older than Hyades stars of the same 
Teff, and this argument is corroborated fe.g. JPace &: Pasauinill2004l) by th e smal l vsini we measure 



(see Table [7|. For TrES-4, the age inferred from the ISestito &: RandichI (|2005l ) relations is t > 5 
Gyr, which is in qualitative agreement with the absence of significant chromospheric activity as well 
as with the inferred value of the surface gravity (see Table [8]) , which indicates the st ar is somewhat 



evolve d. The measured Li abundance for TrES-3 is in line with values determined bv llsraelian et al 



(j2004l ) for a sub-sample of nearby planet hosts with the same Tgfj, while TrES-4 appears more 
depleted than other planet hosts of similar Tefj, \ogg, and [Fe/H]. We discuss briefly in §[7] some 
implications of these findings in the broader context of the existence of chemical peculiarities in 



^One can appreciate how sensitive the age can be to activity cycles by considering the temporal evolution of the 
activity levels in the Sun, with values of log R'h k that range fr om —5.10 to —4.75. These correspond to ages of 
~ 8 Gyr and ~ 2.5 Gyr, respectively (see, e.g.. 



Henrv et alJll99i 
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planet hosts in comparison with stars without detected planets. 

From the measurements in this section we conclude that the above results are consistent with 
the notion that neither star is very young (t > 1-2 Gyr), that TrES-3 is likely to be of intermediate 
age ('^3 Gyr), and that TrES-4 appears older and more evolved. 



4. Revised spectroscopic orbital solution for TrES-3 



Using the radial ve locities presented in ^ 12.11 we have updated the orbital solution given by 
O 'Donovan et al.l (|2007l ). adopting the improved ephemeris described below in ^ 15. li A Kepler ian 
circular orbit was adjusted to the data, a nd the results may be se en in Table [9] and are shown 
graphically in Figure [71 As found also by lO'Donovan et al.l (|2007l ) , the solution gives a scatter 
that is larger than expected from the internal velocity errors. This is most likely due to velocity 
"jitter" associated with chromospheric activity. If we model this as excess scatter to be added 
quadratically to the internal errors, we find that the amount of jitter required to produce a reduced 
value near unity is ~18.5 ms~^. External estimates of the jitter for TrES-3 can be made on 
the basis of the spectral type, the measured value of usini = 1.5 ± 1.0 kms~^, and the activity 
i ndex log-Rirp^ = —4.54 ±0.13, and vary considerably but generally ran ge from about 5 to 20 ms^^ 



(jSaar et al.lll998l : I Santos et al 



200C 



Paulson et al 



2002 



Wrightll2005l ). These estimates are thus 



consistent with our findings. Nevertheless, as a test we also modeled the data with an eccentric 
orbit but did not obtain much improvement, and the result for the eccentricity was not significantly 
different from zero (e = 0.015 it 0.019). There are no indications of long-term variations in the 
observations at hand. 



For TrES-4 we adopt in the following the spectroscopic orbit bv lMandushev et al.l (|2007l ). since 
the radial-velocity material has not changed. For clarity, the orbital solution is reported again in 
Table M 



5. Light-curve analysis 

The analysis of the differentia l photometry f or TrE S-3 and TrES-4 was carried out in essentially 
the same manner as described in Torres et al. (j2008l ). We refer the reader to that paper for the 
details. Given that we have multiple light curves for each system, we first determined the individual 
times of transit as described below in order to improve the ephemeris, and we subsequently refined 
the light-curve parameters. 
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5.1. Transit timings 



To determine the center of each measured transit event we initially adopted the light-curve 
parameters from the disc overy papers. For TrES-3, we adopted a quadratic limb darkening (LD) 
law with coefficients from IClaretl (j2004l ). interpolated to the values T^fi = 5650 K, log 5 = 4.40 dex, 
[M/H] = -0.2 dex, and = 2.0 kms~^ All LD coefficients for the various filters are reported in 
Table O We fitted each light curve individually to solve only for the time of transit center Tc and 
also a linear function of time (two parameters) to describe the OOT fiux (the slope in the OOT fiux 
accounts for systematic errors, including differential extinction). The measured mid-transit times 
are presented in Table [H We then fitted a straight line to the central times of the eight transits 
of the form Tc{E) = Tc{0) + E ■ P and derived the following new ephemeris (which we report in 
Table [9D: = 2454185.9104 ± 0.0001 (HJD), P = 1.30618581 ± 0.00000051 days. With a number 
of degrees of freedom = 6, the resulting reduced x^/^ = 5.87 indicates a rather poor fit. The 
transit timing residuals for TrES-3 are shown in Figure [8j There are four outliers at the 2-3. 5(j 
level. This could be seen as evidence suggesting that the period is not constant. Alternatively, the 
errors might have been underestimated. The data available are not enough to draw any significant 
conclusion on the nature of these variations, but clearly additional observations are warranted. 
We leave for a future study the evaluation o f the dynamical interp retation and signific ance of the 



transi t times, as was done, for example, by ISteffen &: Agoll (120051 ) for TrES-1 and by iDiaz et al 
dioosl) for OGLE-TR-lllb. 



Similarly for TrES-4, we re-fitted the three discovery light-curves using the IClaretl (|2004l ) LD 
coefficients listed in Table [5] appropriate for Tefj = 6200 K, logg = 4.00 dex, [M/H]= +0.1 dex, 
and = 2.0 kms^^. The corresponding times of transit center, derived with the same procedure 
adopted for TrES-3, are presented in Table El Given the short time baseline of these observations, 
we did not derive a new ephemeris for the system, but for the purpose of the analysis presented in 
the next Section we simply adopted the Mandushev et al. (2007) values of P and Tc. 



5.2. Light-curve system parameters 

Next, for both TrES-3 and TrES-4 we locked the transit times and OOT baseline functions 
at the values indicated above and then, under the assumption of a circu lar orbit, we fitted all 



light curves simultaneously using the algorithm by iMandel &: Agoll (j2002l ) to derive the relevant 
quantities radius ratio Rp/R-k (where Rp and i?* are the planetary and stellar radius, respectively), 
inclination i, projected separation a/Ri, (where a is the semi-major axis), and impact parameter b = 
acosi/R^. For TrES-3 we found Rp/R^ = 0.1654 ± 0.0018, i = 81?83±0?12, a/R^ = 5.922 ±0.051, 
and b = 0.840 ± 0.010. For TrES-4, we obtained Rp/R^ = 0.09964 ± 0.00086, i = 82?59 ± 0?40, 
a/R^ = 5.93 lb 0.19, and b = 0.766 ± 0.020. These estimates were derived assuming initial values 
for the stellar and planetary masses of = 0.936 Mq and Mp = 1.920 Mjup for TrES-3, and 
= 1.394 Mq and Mp = 0.923 Mjup for TrES-4, although the results are insensitive to these 
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values. 



Recent studies of transiting exoplanets light curves (e.g., ISouthworthl |2008| . and references 
therein) have highlighted some issues related to how much effect different treatments of limb dark- 
ening can have on the light-curve solutions. In particular, these do not appear to be significantly 
affected by the specific choice of the LD law, while fixing the LD coefficients to their theoretical 
values seems to result in significantly smaller uncertainties in the fitted parameters with respect 
to the case in which the LD coefficients are also adjusted during the fitting procedure. The above 
results were obtained keeping all LD coefficients fixed. We decided to investigate the effect of fitting 
for the LD coefficients in the following way. Given the quality of the data available to us, we can 
only fit for one coefficient (or one combination of the two coefficients) . We chose to fix the quadratic 
coefficient U2 and solve for the linear coefficient ui . We found that allowing complete fr eedom in the 
linear coefficient resulted in "unphysical" solutions (e.g., largely negative coefficients). ISouthworth 
(j2008l ) found that the fitted LD coefficients are usually within 0.1-0.2 of t he theoretical Claret 
values, and that is consistent with our experience (see, e.g.. lWinn et al.ll2007l ). We then decided to 
use an a priori constraint enforcing an agreement of ~ 0.2, modifying the merit function as follows: 



X 



E 



fobs 



fcalc 
J 3 



+ 



Ul 



0.2 



(1) 



where 7?^* is the stellar fiux observed at time j, aj its corresponding error, f9<^^'^ is the model 
value, u is the adjustable linear LD coefficient, and ui its theoretical value (appropriate for each 
band-pass). Furthermore, for TrES-3 all three z-band light curves were required to agree on the LD 
parameter, and the same requirement was set on the two z-band light curves for TrES-4. The results 

81?85 ± 0?16, 
: 82?59±0?40, 



for the TrES-3 system parameters in this case were Rp/R^, = 0.1655 it 0.0020, 



and a/R^, = 5.926 ± 0.056. For TrES-4, we obtained Rp/R^ 
and a/R^, = 5.94 ±0.21. 



0.09921 ± 0.00085, i 



In both cases, the agreement between the system parameters derived keeping the LD coeffi- 
cients fixed and those when the LD coefficients are part of the solution was excellent, with only 
a slight increase i n the estimated uncertainties for the latter case. Indeed, as already noted by 
SouthworthI (j2008l ). fixing the LD coefficients at their theoretically predicted values does not ap- 
pear to significantly bias the results. However, in the interest of providing more conservative error 
estimates, we believe that, for the purpose of the analysis of high-quality light-curves such as the 
ones presented in this paper, a procedure that treats LD coefficients as adjustable parameters is 
preferable. Based on the above considerations, for both TrES-3 and TrES-4 we elected to include 
in Tableland Table [TO] the values of the systems parameters obtained from the light-curve analysis 
in the case in which the line ar LD coeffici e nt wa s allowed to fioat. We note that this approach 
provides a departure from the lTorres et al.l (j2008l ) analysis. 
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6. Stellar and planetary parameters 



The revised spectroscopic determinations of Tefj and [Fe/H], along with the new spectroscopic 
orbital solution for TrES-3 and the light curve fits presented above, allow us to refine the determi- 
nation of the stellar and planetary properties for both systems. To es tablish the properties o f the 
parent stars we rely on stellar evolution models by lYi et al.1 (|200ll ) and lDemarque et alJ (120041 ) . We 
make explicit us e of the constra i nt on the stellar density provided by the light-curve quantity a/i?*, 
as described bv ISozzetti et al.1 (j2007l ). The procedure follows closely that given in the previous 
citation. Briefly, we seek the best match (in a sense) between the measured {Teg, [Fe/H], a/i?*} 
and points on a finely interpolated grid of isochrones spanning a wide range of metallicit i es an d 
ages. Two minor improvements in this procedure, described in more detail bv lTorres et al.l (|2008l ). 
have to do with the weighting of each point sampled along the isochrones according to the distance 
in refj-[Fe/H]-a/i?^ space compared to the observed values, and an additional weighting according 
to the likelihood that the star is in a particular evolutionary state. The latter effect is accounted 
for by multiplying the first weight by the expected number density of stars at each location in the 
H-R diagram, according to an adopted Initial Mass Function (which in this case is simply a power 
law with a Salpeter index). The impact of these weighting factors is generally minor. 

The results for the stellar properties of TrES-3 and TrES-4 are presented in Tableland Table[8j 
The distance estimate to TrES-3 accounts for interstellar extinction (Ay = 3.1 • E{B — V) = 
0.22 lb 0.04 mag), as described in ^ 13.2.41 For both stars the surface gravity inferred from the models 
is typically much more accurate than the spectroscopic gravity determination, showing the power of 
the constraint on the luminosity and size of th e star afforde d by a /R-k. The results for both TrES-3 
and TrES-4 supersede those given recently by iTorres et al.l ()2008l ) because of the new photometry 
and radial velocities contributed here in the former case, and the new light-curve solutions for both 
systems. The evolutionary ages for the two stars are qualitatively consistent with our conclusions 
from § 13.31 but are somewhat more uncertain in the case of TrES-3. Fina lly, using the distance 
estimates in ferred in this work, the UC AC2 proper motio n coni ponents (jZacharias et al.l l2004l ) 
reported bv lO'Donovan et al.l (j2007l ) and iMandushev et al.l (j2007l ). and the mean radial velocity 
values i?VTrES-3 = +9.58 ib 0.73 kms~^ and i?VTrES-4 = —16.40 ± 0.19 kms~^ as measured with 
the CfA Digital Speedometers, we obtain Galactic space motion vectors [U, V, W] = [+27.3, +6.7, 
+33.0] kms-^ and [U, V, W] = [-43.9, -39.1, -6.9] kms^^ for TrES-3 and TrES-4, respectively 
(where U is taken to be positive toward the Galactic anticenter). We collect these results along 
with the other properties derived previously in Table [7| and Table [HI 

The planet parameters follow from the stellar properties and the results of the transit light 
curve and spectroscopic orbits, and are presented in Tableland TablefTOl For TrES-3, our planetary 
mass a nd radius are signif icantly larger compared with the determinations in the discovery paper 
and in [Torres et al.l (j2008l ) . This is due in part to the increased mass and radius for the parent 
star, but also to the larger radius ratio based on the new photometry. 
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7. Summary and discussion 

Our detailed spectroscopic analyses of TrES-3 and TrES-4 have yielded accurate values of the 
atmospheric properties (Tefj and [Fe/H]), which are critical for establishing the fundamental prop- 
erties of the hosts. The accuracy of the temperatures is supported by a number of independent 
checks (low-resolution spectroscopy, line-depth ratios, line profiles, color-temperature calibra- 
tions) that gives us confidence that the inferred stellar properties are reliable. We find that TrES-3 
is a main-sequence G dwarf with a metallicity about 1.5 times lower than the Sun's, not a very 
common occurence among exoplanets hosts. TrES-4 is a somewhat evolved late F star that is near- 
ing the end of its main-sequence phase, and is slightly enhanced in its iron content with respect to 
the solar abundance. 

The agreement we find between age indicators for TrES-3 and TrES-4 based on measurements 
of the Ca II activity levels, the lithium abundances, and rotation, and the evolutionary age inferred 
from the models is fair, although, as discussed in ii l3.3l the reliability of empirical age estimates for 
stars that are not young {t > 1 Gyr) is somewhat questionable. The model estimates themselves 
are not without their problems. Nevertheless, we conclude the stars are 1-3 Gyr old. Neither star 
stands out as peculiar when compared with other planet hosts with similar physical properties. 
In TrES-3 the small vsini value and the fact that we can only place an upper limit on the Li 
abundance are consistent with the notion that planet hosts with T^ff similar to the Sun appear to 
rotate m ore slowly and are more Li-depleted than stars without detected planets. As pointed out 



recently (jGonzalezj |2008| . and references therein), this evidence suggests that a planet-forming disk 
may induce additional rotational braking, leading to enhanced mixing in the stellar envelope, which 
in turn accelerates the destruction o f lithium. TrES-4, on the other hand, does not seem consistent 



with the claim by iGonzalea (|2008l ) that hotter planet hosts with Tcs > 6100 K have higher Li 
abundances, possibly due to self-enrichment processes. It is worth keeping in mind, however, that 
these discrepancies may not be significant given th e large spread in the Li abundance for field 



stars with the temperature and mass of TrES-4 (e.g., iLambert Reddyll2004l ). New investigations 
on these issues are clearly needed based on uniform analyses of large samples of planet hosts and 
statistically significant, well-defined control samples of stars without detected planets. 

New radial- velocity measurements for TrES-3 presented here have enabled us to revise the 
spectroscopic orbit for that system. We detect no indication of any longer-term variations in 
the radial velocities that might suggest the presence of another body in the system. However, 
the small number of observations and their limited time span of only 6 months emphasize the 
need for continued Doppler monitoring of this and other transiting planet systems to investigate 
the possibility of additional companions. In stark contrast to the considerable ground-based and 
space-based efforts invested in studying in great detail the atmospheric properties of many of these 
objects, which have undoubtedly led to tremendous insights into their structure, formation, and 
evolution, the amount of radial- velocity data available for transiting planets is meager, and often 
does not go beyond the handful of observations published in the discovery papers. Interest in the 
radial velocities seems to be quickly lost. It should be pointed out that the frequency of close-in 



- 15 - 



giant planets (P < 10 day) with additional massive planets in outer or bits (up to the detection 



limit of today's Doppler surveys, ~ 4 AU; see, e.g., iButler et al.l |2006| . and references therein) 
is about 12% (8 out of 70 systems discovered via RV methods), and it would therefore be wise 
to continue the velocity monitoring of some of the transiting systems. If additional planets in a 
transiting system were detected, they might also be found to undergo transits, and such a discovery 
would allow us to constrain structural models for gas giants akin to Jupiter or Saturn, and open 
exciting opportunities for additional investigations with present (Spitzer) and upcoming (JWST) 
space-borne observatories. 

Accurate stell ar properties for TrE S-3 and TrES-4 have be en derived here following the ap- 
proach described in lSozzetti et al.l (j2007l ) and lHolman et al.l (|2007l ). comparing the spectroscopically 
determined Tefj and [Fe/H] values along with the photometrically measured a/R^ with current stel- 
lar evolution models. These properties have in turn allowed us to refine the determination of the 
mass and radius of the planets. In particular, we confirm that TrES-4 is the planet with the largest 
radius among the currently known transiting hot Jupiters. 

Recently, Winn et al. (2008a) derived upper limits on the albedo of TrES-3 based on the 
nondetection of occultations observed at optical wavelengths. Our findings are relevant to this 
study in two ways. Firstly, they strengthen the case for a circular orbit, which is important because 
if the orbit is eccentric then it would be possible that Winn et al. (2008a) did not observe TrES-3 
at the actual times of occultations and that their data place no constraint on the albedo. Secondly, 
our revised light-curve parameters are relevant because the upper limit on the geometric albedo 
(pa) was inferred from the measured upper limit on the planet-to-star flux ratio (eA) according to 



Px = ex{a/Rp 



(2) 



Our revised value of (a/Rp) therefore leads to revised upper limits on the geometric albedo of 
TrES-3. However, this revision turns out to be minor: the new value of (a/Rp) is only 1.5% larger 
than the value used by Winn et al. (2008a). The upper limits on the geometric albedo become 
weaker by about 3%. At 99% confidence, the revised upper limits are 0.31, 0.64, and 1.10 in i—, 
z—, and iZ— band, respectively. 

There is a large spread in the observed radii and densities for transiting planets of comparable 
mass placed at similar orbital distance from stars of very similar properties (Tefj, logs', [Fe/H], 
and age). For example, if we consider TrES-4 along with 5 other transiting planet hosts (exclud- 
ing HAT-P-2) with similar characteristics (HD 149026, HD 209458, OGLE-TR-56, OGLE-TR-132, 
and WASP-1), the nominal masses of the attending planets vary by a factor of ^ 3.5, but reported 



densities vary by a factor of ~ 7.5 (e.g., [Torres et al.ll2008l ). Many theoretical mechanisms have 
been proposed to inflate the radius of a strongly irradiated planet, such as additional sources of 



internal heating due to stellar insolation ([Guillot &: Showmanll2002l ). tid al heating due to non-zero 



ml 



eccentricity caused by gra vitational interaction w ith an outer companion (iBodenheimer et al.ll2001 



or by rotational obliquity (jWinn fc Holmanll2005l ). elevated interior opacity due to enhanced atmo- 
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spher ic metallicity (|Guillot et al.ll2006l : iBurrows et al.ll2007l ). or varying core masses (jFortnev et al 
20071 ) . N one of these appear c apable of explaining the observed spread in density in a natural way 
(see also iFabrvckv et alj 120071 ). 



Of the 6 stars just mentioned, all are more metal-ric h than the Sun, except for HP 209458 , 
which has [Fe/H] close to solar. Some of the above models (jGuillot et al.ll2006l : iBurrows et al.ll2007l ) 
predict a positive correlation between the inferred planetary core mass and the host star's metal- 



licity , as in the frainewor k of the core-accretion model of giant planet formation (e.g.. iPollack et al 



19961 : lAlibert et al.ll2005l ). This idea assumes [Fe/H] closely tracks the metallicity of the protoplan- 
etary disk, so that more metal-rich stars should be orbited by more metal-rich planets, with a larger 
heavy-element content. However, among these 6 systems only one planet (HD 149026b) has an in- 
ferred core mass significantly larger than zero, and four of the other planets have measured radii so 
large that the cores are likely to be insignificant. In fact, even in the absence of a core the observed 
radii cannot be reproduced by the models, with TrE S-4 being the extrem e case. Interestingly, over 
40% of the transiting planets reported in Table 5 of lTorres et al.l (120081') do not appe ar to require 
any core at all to explain their sizes (according to the models of iFortney et al.l 120071 ). Taking all 
this into consideration, the claimed evidence for a core mass-met allicity correlation could indeed 
be seen as supporting the more widely accepted scenario of formation by core accretion, but from 
the indica tions above it may a lso be that a significant fraction of these objects formed in a different 
way (e.g., iDurisen et al.l 120071 . and references therein). 



Given the evidence collected so far, we suggest that simply connecting the host star's char- 
acteristics to the structural properties of transiting planets may in fact be an over-simplification. 
We conclude by stressing the importance of refining our understanding of the complex interplay 
between the disk environment and a forming giant planet, and its evolutionary history after en- 
velope accretion, which might turn out to be more directly responsible for its final structure and 
composition than the metal content of the parent star. 
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Table 1. Radial velocity measurements of TrES-3. 





Radial Velocity 


0"RV 


BJD-2,400,000 


(m s~^) 


(m s^"*^) 


54187.04136 


176.0 


8.2 


54187.13832 


99.4 


10.7 


54188.02349 


189.5 


11.2 


54188.12067 


292.4 


13.3 


54189.01696 


-331.6 


5.5 


54189.09680 


-246.4 


11.4 


54189.14038 


-195.8 


7.3 


54288.85394 


241.6 


7.8 


54288.98956 


116.7 


11.5 


54289.82337 


-5.9 


10.8 


54372.86102 


-335.9 


9.5 
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Table 2. Main characteristics of the TrES-3 and TrES-4 hght-curves 



UT Date 


Filter 


Observatory 


Cadence (min) 


OOT RMS 


Res RMS 










TrES-3 








2007-03-25 


z 


FLWO 


1.73 


0.0012 


0.0011 


1.17 


2007-04-08 


B 


FTN 


1.15 


0.0011 


0.0014 


1.46 


2007-04-24 


V 


OT 


1.34 


0.0016 


0.0016 


1.30 


2007-04-25 


9 


FLWO 


0.98 


0.0014 


0.0015 


1.16 


2008-03-09 


i 


FLWO 


0.73 


0.0018 


0.0016 


1.00 


2008-03-27 


r 


FLWO 


0.73 


0.0014 


0.0014 


1.07 


2008-04-12 


i 


FLWO 


1.05 


0.0015 


0.0015 


1.02 


2008-05-08 


i 


FLWO 


0.73 


0.0021 


0.0020 


1.16 








TrES-4 








2007-05-03 


z 


FLWO 


0.73 


0.0015 


0.0016 


1.26 


2007-05-10 


z 


FLWO 


0.73 


0.0019 


0.0017 


1.23 


2007-05-10 


B 


Lowell 


1.54 


0.0015 


0.0015 


1.22 



Note. — Column 4 gives the median spacing between exposures, in minutes. Column 5 
gives the out-of-transit root-mean-squared relative flux. Column 6 gives the residual RMS 
relative flux after subtracting the best-fitting model. Column 7 gives the scaling factor /3 
that was applied to the single-point fiux uncertainties to account for red noise (see ^ 12. 2p . 



Table 3. Differential pfiotometry of TrES-3 



HJD Relative flux Uncertainty 



z band (FLWO 1.2-m) 

2454185.850884 0.99965 0.00115 

2454185.853303 0.99948 0.00115 

2454185.854495 0.99852 0.00115 

B band (FTN 2.0-m) 

2454198.948010 0.99987 0.00109 

2454198.948809 1.00054 0.00109 

2454198.949606 0.99878 0.00109 

V band (OT 0.8-m) 

2454214.574418 1.00159 0.00157 

2454214.575348 1.00247 0.00157 

2454214.576280 0.99766 0.00157 

g band (FLWO 1.2-m) 

2454215.850393 0.99951 0.00150 

2454215.851076 0.99912 0.00150 

2454215.851771 0.99842 0.00150 

i band (FLWO 1.2-m) 

2454535.897759 0.99906 0.00163 

2454535.898268 0.99996 0.00163 

2454535.898766 0.99755 0.00163 

r band (FLWO 1.2-m) 

2454552.866642 0.99987 0.00135 

2454552.867140 1.00181 0.00135 

2454552.867638 1.00086 0.00135 

i band (FLWO 1.2-m) 

2454569.875280 1.00116 0.00149 

2454569.875778 0.99975 0.00149 

2454569.876287 1.00090 0.00149 

i band (FLWO 1.2-m) 

2454594.713269 1.00015 0.00204 

2454594.713778 0.99628 0.00204 

2454594.714276 1.00394 0.00204 
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Note. — The time stamps represent the Helio- 
centric JuUan Date at the time of mid-exposure. 
The data have been corrected for residual ex- 
tinction effects, and the uncertainties have been 
rescaled as described in § [2l We intend for this 
table to appear in its entirety in the electronic 
version of the journal. A portion is shown here to 
illustrate its format. The data are also available 
from the authors upon request. 
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Table 4. Differential photometry of TrES-4 



HJD Relative flux Uncertainty 



z band (FLWO 1.2-m) 



2454223, 


,741853 


0.99446 


0. 


,00160 


2454223, 


,745823 


0.99218 


0, 


,00160 


2454223, 


,746321 


0.99209 


0, 


,00160 




z band 


(FLWO 1.2-m) 






2454230, 


,705029 


0.99892 


0, 


,00170 


2454230, 


,705550 


1.00057 


0, 


,00170 


2454230, 


,706059 


0.99974 


0, 


,00170 




B band 


(Lowell 0.8-m) 






2454230, 


,770620 


0.99944 


0, 


,00150 


2454230, 


,773240 


0.99791 


0, 


,00150 


2454230, 


,774280 


0.99852 


0, 


,00150 



Note. — The time stamps represent the Helio- 
centric Julian Date at the time of mid-exposure. 
The data have been corrected for residual ex- 
tinction effects, and the uncertainties have been 
rescaled as described in § [2l We intend for this 
table to appear in its entirety in the electronic 
version of the journal. A portion is shown here to 
illustrate its format. The data are also available 
from the authors upon request. 



Table 5. Quadratic LD coefficients adopted for TrES-3 and TrES-4 



Filter Linear coefficient ui Quadratic coefficient U2 



TVES-3 

B 0.6379 0.1792 

V 0.4378 0.2933 

g 0.5535 0.2351 

r 0.3643 0.3178 

i 0.2777 0.3191 

z 0.2179 0.3162 

TrES-4 

z 0.1483 0.3600 

B 0.5377 0.2579 



Note. — The assumed limb-darkening law was In/Io = 
1 - ui{l - n) - U2{1 - . 
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Table 6. Mid-transit times of TrES-3 and TrES-4 



HJD Uncertainty (days) Epoch E 



TrES-3 



2454185, 


.910430 


0.000198 





2454198, 


.973147 


0.000223 


10 


2454214, 


.646298 


0.000280 


22 


2454215, 


.952080 


0.000214 


23 


2454535, 


.968246 


0.000166 


268 


2454552, 


.948971 


0.000147 


281 


2454569, 


.929089 


0.000153 


294 


2454594, 


.745943 


0.000253 


313 






TrES-4 




2454223, 


.797215 


0.000847 





2454230, 


.904913 


0.000656 


7 


2454230, 


.905624 


0.001106 


7 



-29- 



Table 7. Properties of the TrES-3 parent star 



Parameter Value 

Teff (K)^ 

logy'' 

log 5^ 

V sini (km s"-*")^ . . . 

(km s~^)^ 

[Fe/Rf 

(log^W)" 

loge(Li)^ 

(g cm~^)^ 

(Mo)^ 

-R* (Rq)^ 

Age (Gyr)b 

L. (i:©)" 

My (mag)^ 

Distance (pc)^ 

U, V, W (km s-i)^ 



^Determined spectroscopically. 

^Inferred from stellar evolution models using ob- 
servational constraints (see text). 

'^Derived observationally. 

Note. — The value adopted for the solar abun- 
dance of iron is log(A?^Fe/-^H)0 = 7.52 



5650 ± 75 

4.4 ± 0.1 

^•ODO_0.014 

1.5 ± 1.0 

0.85 ± 0.1 
-0.19 ± 0.08 
-4.54 ±0.13 

< 1.0 
2.304 ± 0.066 

829+°-°^^ 
^•-'-o.s 

fi25+°-°^6 
U.DZO_Q 

5.39 ±0.11 
228 ± 12 
[+27.3, +6.7, +33.0] 
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Table 8. Properties of the TrES-4 parent star 



Parameter Value 



Teff (K)^ 


6200 ± 75 


logy'' 


4.0 ± 0.1 


log 5^ 




V sini (km s"-*")^ 


8.5 ± 0.5 


(km s~^)^ 


1.50 ± 0.05 


[Fe/Rf 


+0.14 ± 0.09 


(log^W)" 


-5.11 ±0.15 


loge(Li)^ 


< 1.5 


(g cm~^)^ 


Q 014+0.034 


(Mo)^ 




-R* (Rq)^ 


1 O/|«+0.096 


Age (Gyr)b 


^•^-0.4 


L. (i:©)" 




My (mag)^ 


3.13 ±0.17 


Distance (pc)^ 


492 ± 39 


U, V, W (km s-i)^ . . . . 


[-43.9, -39.1, - 



^Determined spectroscopically. 

^Inferred from stellar evolution models using ob- 
servational constraints (see text). 

'^Derived observationally. 

Note. — The value adopted for the solar abun- 
dance of iron is log(A?^Fe/-^H)0 = 7.52 
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Table 9. Revised spectroscopic orbit and light curve solution for TrES-3, and inferred planet 

parameters 



Parameter Value 



Light curve parameters 

P (days) 1.30618581 (fixed) 

Tc (HJD) 2,454,185.9104 (fixed) 

a/R^ 5.926 ±0.056 

Rp/R^ 0.1655 ± 0.0020 

b = acosi/R^ 0.840 ± 0.010 

i{deg) 81.85 ±0.16 

Spectroscopic parameters 

K {ms-^) 369 ±11 

THiREsMms-^) +369.8 ±7.1 

e (fixed) 

{Mp sin i)/{M^ + Mpfl^ [Mq) . . . 0.001893 ± 0.000058 

rms (ms~^) 22.0 

Planet parameters 

Mp (Mj„p) l-SlOlliilis'^ 

Rp (-Rjup ) ^■^^^—O.QSI 

Pp(gcm-3) 0.994lg:°?| 

a (AU) 0.02282l°:S 

log^p (cgs) 3.425 ±0.019 



^7hires is the center-of-mass velocity for the Keck relative 
velocities 

'^The equatorial radius of Jupiter at 1 bar is -Rjup = 71,492 
km. 
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Table 10. Spectroscopic orbit and revised light curve solution for TrES-4, and inferred planet 

parameters 



Parameter Value 



Light curve parameters 

P (days) 3.553945 (fixed) 

Tc (HJD) 2,454,230.9053 (fixed) 

a/R^ 5.94 ±0.21 

Rp/R^ 0.09921 ± 0.00085 

b = acosi/R^ 0.766 ± 0.020 

i (deg) 82.59 ± 0.40 

Spectroscopic parameters 

K {ms-^) 97.4 ±7.2 

THiREsMms-i) +23.7 ±5.8 

e (fixed) 

{Mp sin i)/{M^ + Mpfl^ [Mq) . . . 0.000698 ± 0.000052 

rms (ms~^) 11.1 

Planet parameters 

Mp (Mj„p) 0.925l°:°^i 

Rp (-Rjup*^ ) l-783lQ;Q8g 

Pv 0.202ig:°| 

a im 0.051051°:^ 

log^p (cgs) 2.858 ±0.046 



^7hires is the center-of-mass velocity for the Keck relative 
velocities 

'^The equatorial radius of Jupiter at 1 bar is -Rjup = 71,492 
km. 
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Fig. 1. — Relative flux of the TrES-3 system as a function of time from the center of transit, 
adopting the ephemeris in Table [H Each of the light-curves is labeled with the telescope and 
filter employed. We have overplotted the simultaneous best-fit solution, adopting the appropriate 
quadratic limb-darkening parameters for each band pass (see text for details). 
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Fig. 2. — Relative flux of the TrES-4 system as a function of time from the center of transit, 
adopting the ephemeris in Table [101 Each of the light-curves is labeled with the telescope and 
filter employed. As in Figure [H we have overplotted the simultaneous best-fit model, adopting the 
appropriate quadratic limb-darkening parameters for each band pass (see text for details). 




Fig. 3. — Top: EW of selected Fe I lines in the TrES-3 template spectrum measured manually with 
IRAF vs. EWs measured automatically with ARES. Center: the fractional difference between the 
two measurements as a function of EW. Bottom: Histogram of the fractional differences. 
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Fig. 4. — Observed Hq, profile in the Keck template spectrum of TrES-3 (left) and TrES-4 (right) 
compared with four synthetic spectra with [m/H] = 0.0, log 5 = 4.5, and effective temperatures of 
5500, 5750, 6000, and 6250 K, respectively. 
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Wavelength (Angstrom) 

Fig. 5.— Top: A 20 A region of the Keck template spectrum of TrES-3 centered on the Ca II H 
hne. Bottom: The same, but for TrES-4. 
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TrES — 3 Spectrum 
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TrES — 4 Spectrum 
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Fig. 6. — Top: A 10 A region of the Keck template spectrum of TrES-3 containing the Li I hne at 
6707.8 A (filled dots), compared to three synthetic profiles (lines of various colors and styles), each 
differing only in the lithium abundance assumed. Bottom: The same, but for TrES-4. 
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Fig. 7. — Revised spectroscopic orbital solution for TrES-3, with the post-fit residuals shown at 
the bottom. 
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Fig. 8. — Timing residuals (observed - calculated) for eight observed transits of TrES-3, according 
to the ephemeris derived in this work. 



